SFB/CPP-04-72 
TTP04-28 
hep-ph/0412350 

Strange Quark Mass from Tau Lepton Decays with 0(af) Accuracy 

P. A. Baikov 

Institute of Nuclear Physics, Moscow State University, Moscow 119899, Russia 

K. G. Chetyrkin and J. H. Kiihn 

Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 



in 
o 
o 

(N 

>v 
a ■ 

5 

00 



(N 
> 
O 

in 

m 

CN 

o , 

a: 
i 

^ ■ 



X 



The first complete calculation of the quadratic quark mass correction to the correlator of the two 
currents relevant for the strangeness-changing semihadronic tau-decay rate is presented including 
its real part at the four-loop level. This allows to perform the extraction of the strange quark mass 
m s from the decay width of the tau-lepton with full 0(o? s ) accuracy. In agreement with previous 
estimates the newly computed term proves to be rather large. This justifies inclusion of the 
similarly estimated o4 term in phenomenological analysis. Combined with an updated value of V us 
and an '"improved" version of the renormalization group (RG) improvement of the perturbative 
series this leads to an increase of the central value of m s by about 20% and a partial reduction of 
the theoretical uncertainty by about 50%. 



PACS numbers: 12.38.-t 12.38.Bx 13.65+i 12.20.m 



INTRODUCTION 

The investigation of Cabbibo suppressed semilcptonic 
t decays has developed into one of the important themes 
of r-lepton physics. Combined theoretical and experi- 
mental studies provide an independent and fairly pre- 
cise value for the strange quark mass m s and may 
in the future also lead to a competitive determination 
of the Cabbibo-Kabayashi-Maskawa matrix element V us 
[HESSHIlSlili- The decay proceeds into vector 
and axial current induced final states, which can be fur- 
ther separated experimentally into the spin zero and spin 
one components [lOj. In addition to the total rate also 
various moments of the distribution in s, the invariant 
mass of the hadronic state, can be considered and will be 
discussed below. 

The theory prediction for the rate and moments is 
based on perturbative QCD with small contributions 
from nonperturbative condensates The calculation 
is greatly simplified by the smallness of the strange quark 
mass which justifies an expansion in powers of m 2 s /M 2 . 
Up to now the m s dependence of the total rate was only 
known up to order a 2 . To push the precision further, 
estimates for the third order coefficient were used which 
were based on the principle of minimal sensitivity [l2^ 
(PMS) or fastest apparent convergence [H (FAC). In 
the present paper the exact result for this coefficient is 
presented. It is deduced from the first complete calcula- 
tion of the finite part of a specific correlator in four loop 
approximation. This calculation is based on the concep- 
tual developments described in and requires 
extensive use of computer algebra [17| . The result con- 



firms the PMS/FAC estimate and justifies the use the 
same approach for an estimate of the coefficient. In- 
clusion of this next term leads to a decrease of the central 
value of m s by a few MeV and a partial reduction of the 
theoretical uncertainty by 50%. 

The phenomenological analysis will be based on the, 
by now standard, approach of contour improved pertur- 
bation theory (CIPT) [3 GjJ. Two variants will be dis- 
cussed. We will argue in favour of one of them which 
leads in total to more stable results. The newest devel- 
opments for V us will be included in the analysis. 

The main outcome of the study of the total hadronic 
r width (normalized to the known leptonic width) 
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has been a remarkable confirmation of perturbative 
QCD. It was found that the value of the strong cou- 
pling constant a s as obtained from R T is in good agree- 
ment with those obtained from completely different ex- 
periments such as the Z boson decay into hadrons 
|20l EH E3 . E^ . Furthermore, the strangeness chang- 
ing (Cabbibo-suppressed) part Rs of the decomposition 
of the decay rate and the moments 0] 
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into non-strange (NS) and strange (S) components can 
be used to determine the strange quark mass, one of the 
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fundamental parameters of the Standard Model. The 
moments as introduced above can be experimentally de- 
termined from the measured distribution in the invariant 
mass of the final state hadrons. Hadronic physics is en- 
coded in the quantities r\\ defined through 



of power-suppressed contributions, see |29j). The small 
non-perturbative terms, as well as the m\ terms will be 
effectively included in the phenomenological analysis at 
the end. As a result one has a convenient decomposition 



Rt,ns + Rt.s — 3 Sew (\V u d\ 
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where stands for the CKM matrix and Sew f° r t ne 
universal electroweak correction |2J, |25| . The functions 
r k j are directly related to the correlator of the charged 
current j^x) = ^7^(1 - 7g)V>j 

i f dze^<T[j M (aOit(o)]) = 



(-g^q 2 + g / w t ,)ngV) + 
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through 



k,l = (q)k,l ± r (0)k,X 
ij ij 

2in <f> 
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ds 
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where x = s/M 2 , and the weight functions w^ q ' k ' 1 — 
(1 + 2x)(l - x) k + 2 x l , w^ k ' 1 = -2x l+1 (l - x) k + 2 . The 
behaviour of the perturbative series for the part of the in- 
tegral arising from is more stable than the one from 
However, the latter can either be modelled the- 
oretically and phenomenologically on the basis of scalar 
and pseudoscalar resonance physics 0, H^l i or i be- 

ing solely determined by spin zero contribution in the r 
decays, determined experimentally through the analysis 
of angular distribution of the Kit and K-Kjt, thus separat- 
ing spin zero and spin one contributions [lOj . Restricting 
to scalar and pseudoscalar channels one finds 
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A specific weighted integral over scalar spectral function 
thus directly determines the moments rj° . (A closely 
related discussion along these lines, which deals with 
,/ = 1 and J = spectral functions and the resulting mo- 
ments separately and the corresponding predictions can 
be found in Q.) Since m u ,md <C m 8 and m s <C M T , the 
function r™y is well described by setting m u = m,d = 
and keeping only the leading and quadratic contributions 
to the correlator J3J, viz. 
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Here Q 2 = — q 2 . In the following, we limit ourselves by 
the perturbative contributions (for a recent discussion 
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Let us in addition define the difference 
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which is a useful combination to probe the SU (3) break- 



ing effects as 5; 
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vanishes in the limit 



of exact SU(3) flavour symmetry. In the spirit of the 
decomposition JSJ one similarly defines related quanti- 
ties SR ( T q}kl , 6R { ° ]kl , Sr^ kL and Sr? )kl such that SRf = 
SR ( T q)kl + SR? )kl , Sr kl = Sr^ kl + 5r^ kl . 

The currently known fixed order perturbative predic- 
tions for rn, and 6 u d can be shortly summarized as 
follows Ull EI: 
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(13) 



where a s = a s (M T )/7r. Now, with a s « .1 one observes 
that the "apparent convergency" of the series is accept- 
able for ro but should be considered at best as marginal 
for S us . Clearly, the next term in (jl2(l it is important for 
the interpretation of the measurements. Partial results 
for the term in eq. i|ll|) and the term in eq. (|12fl have 
been published |3l|. We give the results of the calcula- 
tion of the missing term as well as its phenomenological 
implications. 



CALCULATION AND RESULTS 



To compute the real part of {q 2 ) we proceed as fol- 
lows. First, u sing the criterion of irreducibility of Feyn- 
man integrals [16| , the set of irreducible integrals involved 
in the problem was constructed. Second, the coefficients 
multiplying these integrals were calculated as series in 
the l/D — > expansion. Third, the exact answer, i.e. a 
rational function of D, was reconstructed from this ex- 
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TABLE I: Contributions of successive orders of in a s to 
8r lq)kl for the value of a s (M T ) = .334, and normalized to 
the value of 8r^ kl in the Born approximation. First five 
lines: fixed order perturbation theory , second five lines: con- 
tour improved version, with RG improvement for the Adler 
function D 1 - 9 ' ; last five lines: contour improved version with 
RG summation made for the polarization operator IT 9 '. 



(kl) 


Perturbative series 






(0,0) 


1 + 0.425 + 0.283 + 0.178 + 


0.0987 = 


1.98 


(1,0) 


1 + 0.532 + 0.458 + 0.423 + 


0.437 = 


2.85 




i -i- n finfi -i- n ^sq 4- n ro^ 

i u.uuu u.ooy u.u^o ^ 




o.oo 


(3,0) 


1 + 0.663 + 0.695 + 0.793 + 


1.00 = 


4.15 


(4,0) 


1 + 0.708 + 0.784 + 0.943 + 


1.24 = 


4.68 


(0,0) 


0.753 + 0.214 + 0.065 - 0.0611 


- 0.213 


= 0.76 


(1,0) 


0.912 + 0.334 + 0.192 + 0.0675 - 


- 0.0969 


= 1.41 


(2,0) 


1.05 + 0.451 + 0.33 + 0.228 + 


0.0802 = 


= 2.14 


(3,0) 


1.19 + 0.571 + 0.484 + 0.425 


+ 0.33 


= 3. 


(4,0) 


1.32 + 0.697 + 0.657 + 0.665 -4 


- 0.664 = 


= 4.01 


(0,0) 


0.857 + 0.122 - 0.0090 - 0.158 - 


- 0.355 


= 0.458 


(1,0) 


1.11 + 0.232 + 0.11 - 0.0417 - 


- 0.280 = 


= 1.13 


(2,0) 


1.35 + 0.347 + 0.251 + 0.124 - 


- 0.121 = 


= 1.95 


(3,0) 


1.59 + 0.471 + 0.42 + 0.35 + 


0.145 = 


2.97 


(4,0) 


1.83 + 0.61 + 0.623 + 0.648 + 


0.544 = 


4.26 



pansion. Our results read: 
28 



ng s = -4- y a s +a; 
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3 / 6392 



32 
~3 

4496 



(15) 



C3-I6C; 



320 
~27 



27 27 
-4 (-0.35 - 6.437 af) . 

PHENOMENOLOGY 



First of all, it is instructive to compare the exact re- 
sult for the 0(al) contribution to (|14fl with the recently 
obtained predictions ( |3l| ) based on the optimization 
schemes as PMS and FAC: 



d?)3 
v 2,us 



202.309 (exact), 201 (PMS), 199 (FAC). (16) 



This astonishingly good agreement (a similar phe- 
nomenon has been observed for the prediction of the 
Q(a%) term for the correlator of the diagonal currents 
[32J can be considered as a strong argument to repeat 



the procedure and predict, starting from the now com- 
pletely known k^ 3 , the corresponding result for one loop 
more, that is for k^ 4 . To be definite, we use the PMS 



predictions (again for n/ 
similar) 



3; the FAC result is very 



>»4 
v 2,us 



2276 ±200 and k 



(9)4 Aq)i 



2,us 



v 2,itd 



2378±200. (17) 



It is, of course, difficult to estimate uncertainty in the 
above predictions; however, the simple comparison with 
cq. <|16l) clearly demonstrates that an error of about 10% 
should be considered quite conservative. 

Apart of evaluations using fixed order perturbation 
theory two formally equivalent versions of the contour im- 
proved procedure can be found in the literature. The first 
[5( is based directly on the integration of the polarization 
function n^, the second is based on the integration 
of the Adler function = s-j-TL^ and is obtained 

from the first one by partial integration. Partial inte- 
gration and renormalization group improvement do not 
commute as long as finite orders are considered. Since 
the second procedure moves part of the lower order in- 
put to higher orders (contrary to the spirit of CIPT) and 
since, furthermore, the first procedure leads to a some- 
what more stable perturbation series (35|, we consider 
the first of the two choices as preferable. The three op- 
tions for the perturbative series are displayed in Table [i] 
As a consequence of the large value of a s and the rapidly 
growing coefficients of the perturbative series it seems at 
first glance difficult to consider any of the theoretical pre- 
dictions as truly preferable. Nevertheless, the results for 
moments (2,0), (3,0) and (4,0) are at least in plausible 
agreement among the three methods and exhibit accept- 
ably decreasing subsequent terms. Since these moments 
are also relatively most precise, as far as experiment is 
concerned, they will be used in the subsequent analysis. 
The phenomenological analysis will be based on the most 
recent evaluation |33| of \V US \. For the phenomenologi- 
cal description of the contribution due to A°^ kl we adopt 
the analysis presented in [2^ for the second version of 
contour improvement. 

The results for m s , derived from different moments 
and different ways of implementing the contour improve- 
ment procedure are shown in Table [H] (Details about 
the error estimates and the corresponding analysis will 
be given elsewhere.) The main difference, compared to 
the previous analysis, is a downward shift of m s by about 
5 MeV from the inclusion of the aj. terms and an upward 
shift by as much as 20 MeV from the new input for | V us \ . 
The ambiguity for the determined value of m s , including 
the newly computed a\ term, and the estimate for the 
q4 term is shown in Table [H] The renormalzation scale 
H = £M T is allowed to vary between £ = 1 — 1.5. Values of 
fi lower than M T lead to a blow up of a s and destabilize 
the result. By "others" we mean all uncertainties (added 
in quadrature) of the input parameters different from the 
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0(aj?) (or 0(af)) terms in the perturbative contribution. 
They include experimental errors in the moments as re- 
ported in 0, in \V US \ = 0.2259(23) from J33f as well as 
uncertainties in the parameters related to construction of 
the subtracted longitudinal part. We have computed the 
uncertainties using numbers from |28| . 



Parameter 


Value 


(2,0) 


(3,0) 


(4,0) 


m B (0(a^), exact) 




123. 


103. 


88. 


0{a 3 3 ) 


2xO(aJ) 


-3.5 
3.9 


-5.8 
7. 


-6.8 
8.9 


t 


1.5 

1 


-1.4 




4.5 



7.5 



a 3 (M T ) 


.334 ± 0.022 


3.8 
-1.4 


0.44 
1.2 


-1.5 
2.8 


others [2, 28, 33] 




+22.3 
-25.7 


+ 17. 
-19.1 


+ 14.3 
-15.8 


Total 




+23. 
-26. 


+ 19. 

-20. 


+ 18.6 

-17.3 


m 3 (C(4),PMS) 




127. 


100. 


82.4 




2xO(aJ) 
no O(af) 


3.9 
-3.6 


-2.3 
2.4 


-4.6 
5.6 


c 


1.5 

1 


-9.1 




-2. 




4.4 



a a (M T ) 


.334 ± 0.022 


13. 
-5.6 


4.3 
-0.71 


0.24 
2.1 


others [2, 28, 33] 




+22.9 
-26.4 


+ 16.7 

-18.8 


+ 13.8 

-15.3 


Total 




+26.7 
-28.6 


+ 17.5 
-19. 


+ 15.7 
-16. 



TABLE II: Result for m s , derived from different levels of 
approximation, based on contour improvement from [3~H and 
a list of different contributions to the associated error. 

From Tabic □ we find that the inclusion of the a* term 
leads to a better agreement between predictions based 
on two different methods of implementing of the "con- 
tour improvement" approach for the third and the fourth 
moments. As the former moment also shows smaller the- 
oretical error involved we choose it to derive our final 
result for m s which will be given below. 

In total we find 

m s (M T ) = 100+(ii;) theo + (t&) rest MeV . (18) 

If one compares (|18l) to the O(al) result of 28] (cor- 
rected for a different \V US \) m s (M T ) = 106 MeV (with 
larger theoretical and identical remaining errors) one sees 
an essential but still not too large sensitivity to the 
contribution. In fact, for the third moment the shift in 
m s due inclusion of the at* term (-2.5 MeV) is about a 
third of of the corresponding change (-7 MeV) due to the 
a? s contribution. Thus, the purely theoretical uncertainty 
from not yet computed higher orders could be estimated 
as about 3 MeV. Unfortunately, one can hardly hope that 
the error from not yet computed higher orders in a s could 
be reduced further by means of a direct calculation in any 
foreseeable future. 
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